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^ • Abstract. The two-dimensional measurement of the angular correlation of the positron 

annihilation radiation (2D- AGAR) is a powerful tool to investigate the electronic structure of 

C^ , materials. Here we report on the first results obtained with the new 2D- AGAR spectrometer at 

the Technische Universitat Miinchen (TUM). To get experience in processing and interpreting 
2D-AGAR data, first measurements were made on copper. The obtained data are treated 
' O ■ with standard procedures and compared to theoretical calculations. It is shown that the 

C ' measurements are in good agreement with the calculations and that the Fermi surface can 

^ , be entirely reconstructed using three projections only. 



C^ 



1. Introduction 



> 

^fv ' When a positron annihilates in a solid with an electron the resulting annihilation radiation carries 

lO ■ the momentum information of the electron-positron pair. This information is contained in a 

T^lj- I Doppler shift from the mean value of 511 keV, and the angular deviation from exact antiparallel 

^D ■ directions of the two 7 quanta. With 2D-ACAR the latter is measured, which yields a specific 

^^ , projection of the two 7 momentum density p'^'* . Since a positron thermalizes within a few ps [1] 

the angular deviation is caused mainly by the electron momentum in the solid. 

The low momenta that are measured with 2D-ACAR are due to occupied states near the 
^ I Fermi surface in the case of a metal. In these states, any electron is in a delocalised Bloch state. 

J^ ' If we represent the according wave function in the basis of momentum eigenfunctions, it has to 

be expressed by a sum over different momentum components. Mathematically p'^^ is given by 



p27(p) ex V 



i,k:£;(j,k)<£p 



yg-ipr^ep^(r,r)dr 



(1) 



where the sum goes over all the states k below the Fermi level in the jth band. The wave 
function of the positron and an electron in the band j and state k is given by 



*;P,(r,r') = M/^V(r)^^(r')V7.>k(r,r') (2) 



with a general enhancement factor 7j,k(r, r'), which accounts for the interaction of the positron 
with the electron. If we assume an independent particle model, i.e. 7j^k(r, r') = 1, and put in 
the Fourier expansion ^^j^(r)^P(r) = ^q Cj^g exp (i(k + G)r) with reciprocal lattice vectors 
G we get: 

/.27(p) oc Y, 0(^F - i?j,k) E I^G(k) f '^p-k.G (3) 

From this equation we see that filled bands give a continuous distribution and bands crossing 
the Fermi level yield breaks which are distributed through the reciprocal space. An even more 
descriptive result is achieved if we define the following function n(P) of a vector in reciprocal 
space P, assuming a constant positron wave function: 

n(P) = Y.p'''{P + G) oc j;e(i^F - i?,,k) . (4) 

G j,k 

The operation described by equation Q is called the Lock-Crisp- West (LCW) procedure [2]. 
With latter equation the p'^'^ of equation ([3]) can be represented in the reduced zone scheme. 
The LCW procedure is applicable for p'^'^ and its projections in the main symmetry directions, 
i.e. the measured 2D-ACAR spectra. 

The method used here to calculate the full three dimensional p'^'^ is based on the maximum 
entropy algorithm (MEA) [3]. Although a MEA to de-convolve ACAR spectra is known to have 
positive effects on reconstructing p^'^ [4J, it was just recently, that a MEA has successfully been 
applied for reconstructing the 3D momentum density from 2D- ACAR measurements [2- The 
idea behind the MEA can be found, e.g, in |6]. The main purpose of the reconstruction is to get 
a three dimensional density D which agrees within the 2D-ACAR spectra M^ (where i = 1...N'^ 
runs over all pixel) for all measured angles a with the statistical error Uj. It could be realized 
by minimizing the following expression, where Rf{D) is the Radon transformation, i.e. the 
projection, of the density D at an angle of a: 

a i ^ 

If M is of size of N'^, then D has a size of A^'^. Consequently, the minimization problem is 
highly under-determined. This is the reason why the MEA is employed. In order to solve the 
minimization an additional criterion is used: The resulting density D should be flat. For this 
purpose, a Lagrange function with x^ and an entropy function ^ ■ Dj In (Dj) is constructed by 
the use of the Lagrange multiplier A: 

L(Z?) = ^D, In (Z),)- Ax' (6) 

j 
L has to be maximized in order to obtain D. The solution to this problem can be written in an 
iterative form [7] 

with the normalisation constant A. 

In case of p'^'^, the quality of the reconstruction can be improved and the calculation time 
is reduced if the full symmetry of the reciprocal lattice is taken into account. Doing this, the 
MEA becomes a fast and reliable method to recreate the three dimensional p'^^ . 
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(a) Measured (left) and calculated (right) 
anisotropy of the 2D- AGAR spectrum. 



(b) LOW folded density in the [001] 
direction: measurement (left), calculation 
(right bottom) and calculation convolved 
with the experimental resolution (right 
top). 



Figure 1. [001] projection of copper (colour online). 



2. Experiments 

Measurements were performed at the new 2D-ACAR spectrometer at the TUM in Garching. 
The details of the set-up are reported in detail by H. Ceeh et al. [8j. A single crystalline copper 
sample was prepared from commercially available copper single crystal by orienting and cutting. 
To remove lattice defects produced by cutting, the sample was electrolytically polished. 

During the measurement the sample was cooled below 10 K. A spectrum with more than 10® 
counts was recorded for projections along the main symmetry axes [100], [110] and [111]. All 
data were symmetrised according to their point group in order to enhance statistics. For the 
simulation of 2D-ACAR spectra, the Munich SPR-KKR package, version 5.4 by H. Ebert et al. 
was used l9i. 



3. Results and Discussion 



In figure 1(a) the anisotropy of the measured (left) and the calculated data (right) is displayed. 
Obviously the innermost first order features are in very good agreement: There is an increased 
density in [110] direction at the L-point due to the necks of the Fermi surface and a decreased 
density in [100]. The minimum along the latter direction corresponds to the X-point of the 
Brillouin zone. 



This attribute is even more explicit in the LOW folded density on the left side of figure 1(b) 



The highest density is located at the projected L-points, where the Fermi surface extends to 
the boundary of the Brillouin zone, whereas the lowest density is again at the X-points. In the 
case of copper, there are no electronic states at the X-point in the band which crosses the Fermi 



level and hence a low p^'^ is measured. The right side of figure 1(b) displays the theoretical 
LOW density. For a better comparison, the upper right corner shows the calculation convolved 
with the experimental resolution. Again, the congruency is very good. The only small difference 
concerns the size of the necks at the L point. 

To calculate the three dimensional p^'^ a MEA was applied using the three measured 
projections and all the equivalent directions. The data were then folded into the first Brillouin 
zone by the LOW procedure for three dimensions. In figure[2]the isosurface dividing the Brillouin 




Figure 2. Isosurface plot of the reconstructed 3D p^'^ of copper from only three projections. 
The density for the isosurface is chosen in such a way, that the resulting isosurface divides the 
Brillouin zone in two equal parts. 

zone in equal parts is plotted, which has the well known shape of the Fermi surface of copper (see 
e-g- [imfTlJ). Although the qualitative agreement with previous measurements is good, three 
projections are not sufficient to extract reliable quantitative data. 

4. Conclusion and Outlook 

We have performed first measurements with the new 2D-ACAR spectrometer at TUM in 
Garching on the model system copper. The validity of the obtained data could be proven using 
ab-inito calculations together with LOW folding and MEA. Even though only three projections 
were recorded, the 3D p'^'^ could be obtained. 

In the near future, it is planed to measure more projections of copper to compare our 
results quantitatively. Further investigations will then aim at more complex materials, especially 
correlated electron systems. Spectra will be collected as a function of temperature in order to 
investigate phase transitions. In addition we plan to compare the measurements with results 
obtained by complementary techniques such as, e.g., quantum oscillatory methods. 
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